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Abstract 

Similarly to the interaction lagrangian, the possible boundary conditions in quantum field 
theories on space-time manifolds with boundaries are strongly constrained by the symmetry 
and scaling properties of the theory. Based on this general insight, a lattice formulation of 
the QCD Schrodinger functional is proposed for the case where the lattice Dirac operator in 
the bulk of the lattice coincides with the Neuberger-Dirac operator. The construction sat- 
isfies all basic requirements (locality, symmetries, hermiticity) and is suitable for numerical 
simulations. 



1. Introduction 

In numerical lattice QCD, the Schrodinger functional mainly serves as a probe in 
scaling studies of the theory close to the continuum limit [1,2]. The non-perturbative 
computation of renormalization factors, for example, is an important case where the 
Schrodinger functional proved to be very useful (see ref. [3] for an introduction). 

In the continuum limit the Schrodinger functional is obtained by restricting the 
time coordinate Xq to a finite range, < < T, and by imposing Dirichlet boundary 
conditions on the fields at xq = and xq = T. The boundary conditions break chiral 
symmetry, but the chiral Ward identities remain valid away from the boundaries and 
thus forbid the mixing of operator insertions belonging to different chiral multiplets, 
for example, as is the case in the absence of the boundaries. 

The symmetry properties of the Schrodinger functional in lattice QCD should 
ideally be the same if the chosen lattice Dirac operator preserves chiral symmetry 
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on the infinite lattice via the Ginsparg-Wilson relation [4-8]. However, in the pres- 
ence of the boundaries, the Dirac operator and the Ginsparg-Wilson relation must 
both be modified. Chiral symmetry would otherwise not be broken and the lattice 
Schrodinger functional could not possibly have the correct continuum limit. The 
modifications must evidently be local and be linked to the boundary conditions, but 
without further insight it is difficult to say how to proceed from here. 

A theoretically attractive possibility, recently studied by Taniguchi [9] , is to define 
the Schrodinger functional through an orbifold projection. The boundary conditions 
are determined by the Z2 orbifold symmetry in this case, which is taken to be the 
product of a time reflection and a chiral rotation. As it turns out, however, the latter 
leads to some technical difficulties on the lattice and one ends up with a lattice Dirac 
operator whose determinant has a non-removable phase. Moreover, the construction 
becomes rather artificial when the masses of the quarks do not vanish f. 

To a large extent, the solution proposed in this paper is based on universality con- 
siderations. Stated somewhat superficially, the idea is that Schrodinger functional 
boundary conditions do not require any fine-tuning and will thus be satisfied auto- 
matically in the continuum limit, as long as the lattice theory in the presence of the 
boundaries respects locality and the obvious lattice symmetries. The problem then 
reduces to finding a lattice Dirac operator that has these properties and additionally 
satisfies the Ginsparg-Wilson relation away from the lattice boundaries. 

Starting from the Ncubcrgcr Dirac operator in infinite volume [7], such operators 
are not too diflicult to construct. It may be useful, however, to first recall some basic 
facts on the quark sector of the Schrodinger functional (sect. 2) and to address the 
issue of universality in some detail (sect. 3). In sect. 4, Ginsparg-Wilson fermions 
in one dimension are briefly discussed and it is then only a small step to write down 
an acceptable lattice Dirac operator in four dimensions (sect. 5). 



2. Quark fields and Dirac operator in the continuum theory 

In the following, it will be assumed that the reader is familiar with the Schrodinger 
functional in QCD [1,2]. Most of the time the SU(3) gauge field will play a spectator 
role and its presence can be largely ignored. In this section some properties of the 

j" All these problems disappear if chirally rotated boundary conditions are adopted, as suggested by 
Sint [10]. The goal here, however, is to define the Schrodinger functional with the standard parity- 
conserving boundary conditions [2]. 
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quark sector in the continuum theory are highhghted. This serves partly to introduce 
the subject and partly to guide the construction of the lattice theory in the later 
sections. 



2.1 Boundary conditions 

It suffices to consider a single quark with mass m, since all formulae and comments 
trivially extend to the case of several quarks with arbitrary masses. The quark and 
antiquark fields il^{x) and '^{x) are defined at all times xq G [0,T] and are required 
to satisfy X 



P+ip{x) = il){x)P- = at xo = 0, 
P-il){x) = ^{x)P+ =0 at xo=T, 



(2.1) 
(2.2) 



where P± = ^(1 ± 70). These boundary conditions are invariant under space rota- 
tions, parity, time reflections {xq — ^ T — xq) and charge conjugation. 

The quark action in the presence of a gauge field ^^(x) is then given by 



Sf = J dxo J d^xtp{x)D^'ip{x), 
Dm = lixD^j, + m, Di^ = d^ + A^. 



(2.3) 
(2.4) 



Usually the Schrodinger functional is set up with inhomogeneous boundary condi- 
tions, where the boundary values serve as sources for the quark field at a^o = and 

Xq = T [2]. One may, however, just as well adopt homogeneous boundary conditions, 
as is done here, and introduce the boundary quark fields directly through 



ax) = P-^{x), Q{x)=i>{x)P^ at a;o = 0, 
C'(x) = P+V'(x), X' {x) = i^{x)P- at xo = T. 



(2.5) 
(2.6) 



These are just the non-zero Dirac components of the quark fields at the boundaries. 



\ The notation is the one commonly used in lattice QCD. In particular, the space-time metric is 
euclidean and the Dirac matrices 7(_j, /i = 0, . . . , 3, are taken to be hermitian. The fifth Dirac matrix, 
75 = 70717273, is then also hermitian. Where appropriate, the Einstein summation convention is 
applied. 
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2. 2 Spectrum of the Dirac operator 

Let H be the linear space of all smooth quark fields that satisfy the required bound- 
ary conditions. With respect to the natural scalar product in this space, 

(V',x)= dxo J d^xi^ix)'<x{x), (2.7) 

the operator = 75-Dm is easily shown to be hermitian. Since Qm is also elliptic, it 
follows that this operator has a complete orthonormal set of eigenfunctions Vn £ H, 
labelled by an index n = 0, 1, 2, . . ., with real eigenvalues A„. 

For non-negative quark masses m, the eigenvalues are bounded from below by 

A^>m2 + /x2, (2.8) 

where /i denotes the spectral gap at m = 0. This is a straightforward consequence 
of the identity 

= IIQm^^nll^ = rn? \\Vnf + ||QmL=o^n||^ 

+ m d^XVn{x)'^Vn{x) + m / d^XVn{x)^Vn{x). (2.9) 

Jxo=0 J xa=T 

Moreover, using the fact that the current i'n(x)^7^t;„(a;) is conserved at m = A„ = 0, 
it is possible to show that the massless Dirac operator has no zero modes. 

Note that the bound (2.8) becomes invalid at m < 0. The Dirac operator has 
eigenmodes localized at the boundaries in this case, with eigenvalues that decrease 
exponentially at large T. These are just the well-known domain wall fermion modes. 
The fact that the sign of the quark mass matters is a consequence of the breaking of 
chiral symmetry through the boundary conditions. This is also why the eigenvalues 
at non-zero quark masses are not simply related to those at vanishing mass. 

2.3 Determinant of the Dirac operator 

Since there are positive and negative eigenvalues A„, the determinant, dot Qm, may 
have a non-trivial phase. However, this is not the case if the standard lattice regu- 
larization is employed [2]. Other regularizations may give different results, but the 
universality of the continuum limit implies that the determinant must always be 
real after removal of the ultraviolet cutoff, up to some local terms perhaps. These 
can always be "renormalized away" , and the requirement of a real determinant then 
becomes part of the definition of the theory. 
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At non-negative quark masses, the quark determinant actually has a definite sign 
in this case, because the eigenvalues of Qm never pass through zero. In the functional 
integral, the quark determinant may thus be replaced by its absolute value. 

2.4 Chiral symmetry properties of the quark propagator at m = 

As a differential operator, the massless Dirac operator anticommutes with 75 and the 
local chiral Ward idenitites are thus not affected by the presence of the boundaries. 
The quark propagator S{x,y), on the other hand, transforms according to 

75S{x,y) + S{x,y)j5 = 

[ d^zSix,zh,Siz,y)+ [ d^zS{x,z)^,S{z,y). (2.10) 

Jzo=0 Jzo=T 

This equation is easily established by noting that, as a function of x, the expression 

on the left solves the homogenous Dirac equation. The other side of the equation 
is then obtained by reconstructing the solution from its boundary values at = 
and xq = T. 

In QCD with more than one massless quark, eq. (2.10) implies the non-singlet 
chiral Ward identity 

(A'*75V'(^)V^(y) + V'(^)?(y)A"75>F = 

I d^z (V(x)V^(y) {C(;s)A"75C(^) + C'(-2)A"75C'(^)}>p , (2-11) 

where A" is any traceless matrix in flavour space and (. . .)f denotes the expectation 
value in the quark sector. The breaking of chiral symmetry through the boundary 
conditions is made explicit by this formula. Purely from the symmetry point of 
view, it actually looks as if there were a unit mass term at the boundaries. 



3. Field theories on lattices with boundaries 

In continuum field theories, the notion of smoothness plays a central role. Differential 
operators can only act on smooth functions, for example, and boundary conditions 
really only make sense when imposed on such functions. 
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The situation in lattice field theory is quite different. Strictly speaking, boundary 
conditions are no longer imposed on the fields but instead are encoded in the lattice 
action. As usual different lattice theories may have the same continuum limit, where 
now the characterization of the latter must include a specification of the boundary 
conditions. An important point to note is that the possible boundary conditions are 
strongly constrained by the requirement of locality, the lattice symmetries and by 
power-counting arguments. There are then not many more universality classes than 
there are in the absence of boundaries. 

The aim in this section is to make these remarks a bit more concrete, so that they 
can be applied to the Schrodinger functional in QCD. Eventually, the argumentation 
is based on Symanzik's work on the renormalization of quantum field theories in the 
presence of boundaries [11] and also on the theory of boundary critical phenomena 
in statistical mechanics (see ref. [12] for a review). 

3. 1 Boundary conditions and the field equations 

The concepts that will be developed in the following are best introduced by consid- 
ering a free scalar field in the half-space xq > with various boundary conditions 
at xo = 0. A simple lattice formulation of this theory is obtained by choosing the 
lattice field ^(x) to reside on the sites of a hypercubic lattice with spacing a. The 
expression 

= E E i + m'<t>{xf} (3.1) 

xo>a X 

is then a possible choice of the lattice action, where 5^ denotes the forward nearest- 
neighbour difference operator in direction ^ and the mass m is assumed to be pos- 
itive. Note that the action depends on the field variables at time Xq = a, 2a, 3a, . . . 
only. These are thus the unconstrained dynamical degrees of freedom of the field 
which are to be integrated over in the functional integral. 

Starting from the action (3.1), the propagator {(p{x)(f){y)) can easily be worked out 
in a time-momentum representation. It then turns out that the propagator satisfies 
Neumann boundary conditions in the continuum limit. On the other hand, after a 
slight modification of the action by a boundary term, 

S^S + a'^^Jix)\^^^, (3.2) 

X 

the calculation yields a propagator that satisfies Dirichlet boundary conditions in the 
continuum limit, for any fixed c > (the powers of a are such that c is dimensionless). 
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Some understanding of how a particular boundary condition arises can be achieved 
by considering the field equation 



SS 



{r){x)(p{y)) = a "^S^y, r){x) = -jJTZX- (3-3) 



A short calculation yields 

ri{x) = {-d^dn + m'^}(p{x) at xq > a, (3.4) 

where d* denotes the backward nearest-neighbour difference operator. In the bulk of 
the lattice, the field equation thus coincides with the lattice Klein-Gordon equation. 
At xo = a, however, a different equation is obtained, 

c 1 

ri{x) = -^^{x) - -9o(t>{x) + {-dkdk + (3.5) 

which depends on the details of the action close to the boundary. Formally, the first 
term in this expression dominates in the continuum limit and the field equation at 
the boundary thus implies Dirichlet boundary conditions at a = if c > 0. On the 
other hand, the second term dominates if c = 0, which leads to Neumann boundary 
conditions in the continuum limit. 



3.2 Natural boundary conditions 

The example considered in the previous subsection illustrates the fact that bound- 
ary conditions are a property of the continuum theory which arises dynamically in 
the continuum limit. In particular, they arc not determined by the space of lattice 
fields alone. Another outcome is that some boundary conditions are obtained gcner- 
ically, while others (Neumann and mixed boundary conditions in the study case) are 
unstable under perturbations of the lattice action. 

In an interacting theory, and also in free theories with complicated lattice actions, 
the connection between the lattice field equations and the boundary conditions in 
the continuum limit may not be as transparent as suggested above. It seems plau- 
sible, however, that the boundary conditions are always of the form 0{x)\^^_q = 0, 
where 0{x) is a linear combination of local fields with the appropriate number of 
components and symmetry properties. 

The boundary conditions that are stable under perturbations of the lattice theory, 
and which thus arise naturally, then correspond to the fields with the lowest possible 
dimension. Some tuning of the lattice theory will normally be required in all other 
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cases, unless there are some symmetries that protect 0{x) from mixing with lower 
dimensional fields. 

3.3 Are Schrddinger functional boundary conditions natural? 

Following Wilson [13], the quark fields are represented on the lattice by Dirac spinors 
ip{x) that reside on the points x of the lattice. Since QCD is asymptotically free, the 
scaling dimension of the local fields in this theory coincides with their engineering 
dimension (up to logarithms). In particular, the fields of lowest dimension that carry 
the quantum numbers of the quark fields are the quark fields themselves. 

The boundary conditions at xq = 0, which arise naturally in the continuum limit, 
are thus of the form -B^(x)|^^^o = 0) where S is a constant matrix in Dirac and 
colour space. At xo = T the boundary conditions are linked to those at a;o = by 
the time reflection symmetry and charge conjugation then determines the boundary 
conditions for the antiquark field. 

If the lattice theory is invariant under gauge transformations, cubic rotations and 
parity, as will be the case in the following, the boundary conditions must respect 
these symmetries f. Moreover, B cannot have maximal rank, as otherwise the bound- 
ary conditions and the Dirac equation at < Xq < T would imply a vanishing quark 
propagator. Up to a normalization constant, the only matrices that are compatible 
with all these conditions then are B = and B = P^. 

Schrodinger functional boundary conditions thus arise naturally and do not require 
any particular adjustments of the lattice action. There are two classes of lattice theo- 
ries, which arc distinguished by the sign in the boundary condition P±ip{x)\^^^Q = 0. 
The difference matters if the quark mass does not vanish, but the sign may easily be 
determined, in any given case, by studying the free-quark propagator for example. 



4. Ginsparg— Wilson quarks in one dimension 

In the presence of the boundaries, an acceptable lattice Dirac operator that satisfies 
the Ginsparg- Wilson relation in the bulk of the lattice still needs to be found. The 
Wilson-Dirac operator provides a simple solution to this problem in one dimension. 

f Gauge transformations include the gauge field variables at the boundaries and are thus a spurion 
symmetry to some extent. In the present context, this is of no importance, however, since B does 
not depend on the gauge field. 
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This is a somewhat trivial case, but it gives important hints for the construction of 
the Dirac operator in higher dimensions. 

4-1 Lattice Dirac operator 

For simplicity the gauge field is omitted in this section. The Wilson-Dirac operator 
in one dimension then reads 

D = j {70(^0 + do) - ad^do} = P+^o* - P-do. (4.1) 

In the absence of boundaries, this operator satisfies the Ginsparg-Wilson relation 

+ D75 = aDjsD. (4.2) 

Moreover, from the definition (4.1) it is immediate that 75-D is hermitian and that 
the quark propagator is given by 

S{xo,yo) = 0{xo > yo)P+ + 9{xo < yo)P-, (4.3) 

where 9{*) is equal to 1 if the logical condition in brackets is true and otherwise. 

Following the standard treatment [2], the dynamical degrees of freedom of the 
quark fields in the presence of the boundaries at Xq = and Xq = T are taken to be 
their components at a;o = a, 2a, . . . T — a. It is convenient to assume that the fields 
are defined at all other times as well, but that they are equal to zero there. 

The Dirac operator in the presence of the boundaries (which is also denoted by 
D) maps this space of quark fields into itself. At xq = a, 2a, . . . ,T — a, its action is 
given by eq. (4.1), while at xq < and xq > T the target field is set to zero. With 
respect to the scalar product 

T-a 

{ijj, x)=aJ2 V'(a;o)^x(a^o), (4.4) 

xo=a 

the operator 75-D is then again hermitian. It is also easy to show that D has no 
zero modes and that the associated propagator coincides with the propagator (4.3) 
on the infinite lattice, at all times in the range < xo,yo < ^- In particular, the 
correct Schrodinger functional boundary conditions are obtained in the continuum 
limit. 
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4-2 Modified Ginsparg-Wilson relation 

As already pointed out in sect. 1, the lattice Dirac operator is not expected to satisfy 
the Ginsparg-Wilson relation in the presence of the boundaries, for general reasons. 
On the finite lattice, the operator D introduced above in fact satisfies the modified 
relation 

75L> + L>75 = aL>75-D + 75P, (4.5) 
where P is a local operator given by 

PV'(^o) = I {6,,aP-tl^{a) + <5,o T-aP+HT - a)} . (4.6) 

Equation (4.5) is a straightforward consequence of the precise definition of D in the 
presence of the boundaries. In particular, the boundary term 75P is obtained when 
the product -D75-D is worked out at xq = a and xo = T — a. 

By multiplication of the modified Ginsparg-Wilson relation (4.5) from both sides 
with the quark propagator, it follows that 

75'S'(xo, yo) + S{xo, 2/0)75 = 75(^x03/0 + 

S'(xo,a)75P_5(a,yo) + S{xo,T - a)-f5P+S{T - a,yo). (4.7) 

In the continuum limit, the projectors P± on the right-hand side of this equation 
can be dropped in view of the boundary conditions satisfied by the propagator. The 
relation then reduces to the one-dimensional form of the chiral Ward identity (2.10). 
This shows that the presence of the boundary term in the modified Ginsparg-Wilson 
relation (4.5) is directly linked to the breaking of chiral symmetry in the continuum 
theory by the boundary conditions. 



5. Lattice Dirac operator in four dimensions 

Since Schrodinger functional boundary conditions arise naturally, the choice of the 
lattice Dirac operator is not critical and there are probably many viable construc- 
tions. The operator proposed here is a simple modification of the Neuberger-Dirac 
operator in infinite volume. 
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5.1 Definition 

As before the theory is first set up on the infinite lattice. The quark fields are thus 
assumed to be defined at all sites of the lattice. Although the gauge field continues to 
play a spectator role, it is now included in the formulae. The Wilson— Dirac operator 
is then given by [13] 

i^w = H7M(VM + V^)-aVX}. (5-1) 

where V^j and V* denote the gauge-covariant forward and backward difference op- 
erators. As already mentioned, the starting point in this section is the Neuberger- 
Dirac operator [7] 

D = ^{l-A{A^Ay'^^}, (5.2) 

A=l + s-aD^, 0,= -^. (5.3) 

i. ~\~ s 

The parameter s in this formula allows for some optimization and is only relevant in 
the context of numerical simulations. In practice, it is normally set to a fixed value 
in the range < s < 1/2. 

In the presence of the boundaries at = and xq = T, the dynamical degrees of 
freedom of the quark fields reside on the lattice sites at time xq = a,2a, . . . ,T — a. 
It is again convenient to assume that the fields are defined at all other points as well 
and that they are equal to zero there. The Wilson-Dirac operator may be considered 
to be a linear operator in this space of fields, whose action at < xq < T is given by 
eq. (5.1) (elsewhere the target field is set to zero). This is the lattice Dirac operator 
that was introduced by Sint [2]. 

The structure of eq. (5.2) is such that D satisfies the Ginsparg-Wilson relation 
(4.2) automatically (with a replaced by a) if 75^ is hermitian. In the presence of the 
boundaries, the Dirac operator must therefore be given by a different expression. A 
formula that works out is 

D = ^{l-l{U + U~)}, (5.4) 

U = A{A^A + caPy^^^ , U~ = j5U^5, (5.5) 
where c > 1 is another tuneable parameter, whose optimal value will turn out to be 
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close to 1 + s. In eq. (5.5), the boundary operator 



P^(x) = ^ {Sa;oaP-'(p{x)\xo=a + Sxo T-aP+1pix)\xo=T-a} (5.6) 

is the four-dimensional version of the operator P previously encountered, while A is 
again given by eq. (5.3), where is now the Wilson-Dirac operator in the presence 
of the boundaries. 

The merits of the definition (5. 4), (5. 5) will be discussed in detail, but before this 
it may be helpful to note that the operator D reduces to the Wilson-Dirac operator 
in the one-dimensional theory if s = and c = 1. The results reported in sect. 4 
actually imply that A^A + caP = 1 in this case. 

5.2 Lattice symmetries, hermiticity and spectral bounds 

It is not difficult to check that the Dirac operator D transforms like the Wilson-Dirac 
operator under cubic rotations, parity, time-reflections and charge conjugation. The 
latter interchanges U with U~, and having the sum of these two operators in eq. (5.4) 
also ensures that 75 D is hermitian. 

Another implication of the form (5.5) is the bound 

||[/|| = IICZ-II < 1. (5.7) 

To show this, it suffices to note that 

\\U^^ = {x,A'^Ax) < {x,{A'^A + caP)x) = (5.8) 

for any quark field ■0, where x = {A'^A + caP)~^/^ip. The spectrum of aD is thus 
contained in the unit disk in the complex plane centred at 1. However, one should 
not expect the spectrum to be on the unit circle, as is the case for Dirac operators 
satisfying the Ginsparg-Wilson relation. 

In the present framework, the natural choice of the massive Dirac operator is [14] 

Dm = {l-^am)D + m. (5.9) 

When the bare quark mass m is in the range 0<m<2/a, as will be assumed in the 
following, the spectrum of this operator is separated from the origin by a distance 
of at least m. The middle term in the expansion 

{l5Dmf =m^ + m(l - ^am){D^ + D) + {I - \amf {-^^Df (5.10) 
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is in fact non-negative and the eigenvalues A„ of 75-Drre are therefore bounded by 

Xl>m^ + {1- \amfn'^, (5.11) 

where fi denotes the spectral gap at m = 0. This bound coincides with the spectral 
bound (2.8) in the continuum theory, up to corrections of order am. 

An important consequence of these results is that the determinant dct is real 
and non-zero at all quark masses m > 0. Actually, det Dm must be positive at these 
masses, because this is trivially the case at m = 2/ a and because the determinant 
is a continuous function of m, on any finite lattice. 

5.3 Locality 

In position space, the Dirac operator is represented by a kernel D{x,y) through 

T-a 

D^{x) = a^Y.Y. 2^)^(2/), < xo < T. (5.12) 

yo=a y 

Locality requires that the bound 

a^\\D{x,y)\\ < Ce-'^ll^-^ll/" (5.13) 

holds for some constants C and k > that do not depend on a. Moreover, up to such 
exponentially small tails, D{x,y) should be locally constructed and only depend on 
the gauge field variables in the vicinity of x and y. 

In infinite volume, a rigorous proof of the locality of the Neuberger-Dirac operator 
can be given if the gauge field is not too rough on the scale of the lattice spacing [15]. 
Further studies then suggest that locality holds under far more general conditions, 
including those typically encountered in numerical lattice QCD at lattice spacings 
a < 0.1 fm [15,16]. 

The proof presented in ref. [15] is based on an expansion of the inverse square root 
in eq. (5.2) in Legendre polynomials. The expansion converges rapidly if A^A > a 
for some a > 0, and the locality of the Dirac operator then follows immediately. 
In the case of the operator (5.5), the Legendre expansion similarly links its locality 
properties to the existence of a non-zero lower bound on A'^A + caP. 

As shown in appendix A, the spectral gap of A^A+caP is not smaller than that of 
A'^A on the infinite lattice, independently of how precisely the gauge field is extended 
from the time slice < Xq < T to all times. The field may be extended through time 
reflections at the planes xq = mod T, for example, which is a good choice in the 
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present context, since the smoothness properties of the field (if any) are preserved. 

In particular, the estimates of ref. [15] then immediately imply the existence of a 
spectral gap if the gauge field at < < T is sufficiently smooth on the scale of 
the lattice spacing. 

Presumably the locality properties of the Dirac operator (5.4) are thus as good as 
those of the Neuberger-Dirac operator on lattices with periodic boundary conditions, 
also when the rigorous arguments of ref. [15] do not apply. However, some numerical 
studies may still be required to confirm this in the cases of interest. 

5.4 Chiral symmetry 

At a distance d from the boundaries, the kernel D{x, y) of the Dirac operator (5.4) 
coincides with the kernel of the Neuberger-Dirac operator on the infinite lattice, 
up to terms that decrease exponentially like ^-'^'^1°- , From the expansion in Legen- 
dre polynomials mentioned in the previous subsection, for example, this property 
is evident, taking into account the fact that the operator P is supported at the 
boundaries of the lattice. The separation of bulk and boundary terms is actually a 
direct consequence of the locality of the operators involved (see appendix B). 

The remark has two important implications. First of all, it shows that the Schro- 
dinger functional constructed here probes the right theory, i.e. the one where the 
lattice Dirac operator in the absence of the boundaries is equal to the Neuberger- 
Dirac operator. Secondly, it follows that 



where is a local operator with kernel AB(x,y) supported in the vicinity of the 
boundaries (up to the usual exponentially small tails). Correlation functions of local 
fields at physical distances from the boundaries thus satisfy the same chiral Ward 
identities as they do on lattices with periodic boundary conditions, for example. 

Starting from the definition (5. 4), (5. 5) of the Dirac operator, or from the formulae 
in appendix B, the operator can be worked out explicitly. One may hope to find 
Ab = 75-P, as is the case in one dimension, but the expressions that are obtained 
are complicated and not very illuminating. 

5.5 Boundary fields and 0(a) improvement 

A possible lattice representation of the boundary fields (2.5) at time = is 



75L> + L>75 = aD-f^D + Ab, 



(5.14) 



ax) = U{x,0)U,=oP-i^{x)\ 



(5.15) 
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P+U{x,(}) 



-1 



(5.16) 



xo=a 



xo=Oi 



where U {x, fi) denotes the hnk variable at the point x in direction fj,. This definition 

coincides with the one commonly adopted in lattice QCD with Wilson quarks [2]. It 
should be noted, however, that the normalization of these fields depends on the de- 
tails of the lattice regularization and may not be the canonical one (cf. subsect. 5.6). 

Formulations of lattice QCD with Ginsparg- Wilson quarks are automatically on- 
shell 0(a)-improved [17,18]. In the presence of the boundaries, this is no longer the 
case, but the theory can be improved by including a few 0(a) boundary counterterms 
in the lattice action. The list of terms that must be added was determined in ref. [18]. 
One of the counterterms amounts to a modification of the lattice Dirac operator, 
while all others are either pure gauge terms or reduce to contact terms and 1 + (am) 
renormalization factors. 

These other counterterms can be implemented as in the standard lattice theory 
[2,18]. The situation is a bit more tricky in the case of the countcrtcrm that modifies 
the lattice Dirac operator, because some of the desirable properties of the latter may 
be compromised [the spectral bound (5.11), for example]. 

The counterterm at the boundary xq = is usually taken to be a straightforward 
lattice implementation of the boundary action 



The precise choices that one makes do not matter, since the counterterm is uniquely 
determined by its symmetries and dimension, up to redundant terms and corrections 
of higher order in a. In the present context, the 0(a) improvement can therefore 
also be achieved by tuning the coefficient c on which the Dirac operator (5. 4), (5. 5) 
depends. Changes of this coefficient actually amount to a modification of the oper- 
ator in the vicinity of the boundaries by a local term with the correct symmetries. 
The properties of the Dirac operator discussed in the previous subsections are then 
preserved. 

5.6 Free-quark theory 

In the absence of the gauge field, it is possible to check explicitly that the lattice 
theory has the correct continuum limit and that the 0(a) improvement works out 
in the way described in the previous subsection. 




(5.17) 
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The operator under the square root in eq. (5.5) assumes the form 

A^A + caP = (1 + sf + sa^ ^ d*d^ + ^ dld^^- + (c - (5.18) 

At ix<v 

in the free theory. As explained in subsect. 5.1, the operator aets on quark fields in 
the presence of the boundaries. Its eigenfunctions at c = 1, for example, are given 

by 

sin(poXo)e*P"', VQ = —, n = 1, 2, . . . , T/a - 1. (5.19) 

Since the associated eigenvalues are bounded from below by (1 — |s|)^ if \s\ < 1, the 
locality of the free Dirac operator is guaranteed at all |s| < 1 and c > 1. 

In the time-momentum representation, the kernel D(x,y) of the Dirac operator 
can be worked out analytically to some extent. The quark propagator, on the other 
hand, may be difficult to obtain in closed form, but it can be computed numerically 
on lattices with hundreds of points in each direction, using established techniques 
(see ref. [19], for example). This allows the lattice propagator to be compared with 
the continuum propagator in a large range of lattice spacings, spatial momenta and 
quark masses. 

Many checks were then performed, including the following three: 

(1) It was verified that the quark propagator {|^(a;)V'(y)) at non-zero distances from 
the boundary as well as the boundary-to-bulk propagator {C{x)ip(y)'^ have the cor- 
rect continuum limit (explicit expressions for the continuum propagator can be found 
in ref. [20]). 

(2) The full spectrum of (75-D)^ was computed and compared with the spectrum 
in the continuum theory. In particular, the presence of any additional modes that 
would survive in the continuum limit could be excluded in this way. 

(3) The approach of the propagator to the continuum limit was studied and it was 
shown that the lattice effects of order a can be cancelled by tuning the parameter c 
of the lattice Dirac operator and the field normalization factors. 

The 0(a) improvement is achieved for values of c close to I + s, and the improved, 
canonically normalized boundary field is given by Z{1 + ham)C, where Z ~ 1 — s/4 
and h ^ —3/4. These normalization factors will be rarely needed, but it is reassuring 
to note that they are in a range that is not uncommon for such factors. Moreover, 
the observed convergence of the propagators to the continuum limit is quite similar 
to the one seen in the standard Wilson theory [20] . 
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6. Concluding remarks 

Lattice Dirac operators that satisfy the Ginsparg-Wilson relation have little in com- 
mon with the difference operators that are obtained by "discretizing" the classical 
Dirac equation. Classical concepts can in fact be rather misleading in lattice field 
theory. In particular, boundary conditions cannot simply be imposed on the lattice 
fields but arise dynamically in the continuum limit. 

Universality considerations gain additional importance in this context, since they 
apply to any local theory, independently of how complicated it may be. The defi- 
nition of the Schrodinger functional proposed in this paper heavily builds on such 
arguments. It is clearly not the only possible construction, but the proposed formu- 
lation works out and has many attractive features. 

In lattice QCD with Ginsparg-Wilson quarks, the application of domain decom- 
position methods [21] appeared to be excluded so far, because it was not clear how 
to restrict the lattice Dirac operator to blocks of lattice points. Following the lines 
of sect. 5, it is now straightforward to come up with viable expressions for the block 
Dirac operators. The boundary operator P, for example, may be fixed by requiring 
eq. (A.l) to remain valid when R is set to the block restriction operator. The locality 
of the block Dirac operators as well as a number of other desirable properties are 
then guaranteed. 

I am indebted to Stefan Sint for helpful discussions on the Schrodinger functional 
with chirally rotated boundary conditions and to Peter Weisz for a critical reading 
of the paper. 



Appendix A. Spectral bound on A^A + caP 

As explained in subsect. 5.3, the locality of the Dirac operator (5.4) can be proved 
rigorously if A^A + caP > a for some a > 0. In this appendix, it is shown that such 
a bound can be obtained by relating A^A + caP to the operator A^A on the infinite 
lattice. 

The standard setup of the lattice Schrodinger functional assumes the gauge field 
variables U{x, n) to be defined at all times 0<a;o<rif;U = and additionally at 
xq = T ii n = 1,2,3. Evidently, a given field can always be extended to all times by 
setting the so far undefined link variables to unity, for example. 
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Exactly which extension is chosen will not matter in the following. Once a definite 
prescription is adopted, the Wilson-Dirac operator becomes a well-defined, bounded 
linear operator in the space of all square-summable quark fields on the infinite lattice. 
The associated operator A [eq. (5.3)] is denoted by A in order to avoid any confusion 
with the operator A that appears in eq. (5.5). 

A straightforward calculation, similar to the one that leads to eq. (4.5), now shows 
that 

R{A^A + aP) R = RA^AR, (A.l) 

where R projects the quark fields to the physical subspace in the presence of the 
boundaries. Explicitly, 

f ip(x) if < xo < T, 
R^{x) = <^ ^ (A.2) 
[ otherwise, 

for any quark field ip on the infinite lattice. Note that the product on the left-hand 
side of eq. (A.l) is perfectly well-defined, since R projects the fields to the physical 
subspace before the operator in brackets (which can only act on quark fields in this 
space) is applied. 

Quark fields in the physical subspace satisfy Rip = tp, and from eq. (A.l) it then 
follows that 

{iP, {A'' A + caP)V') = {tp, A^AtP) + {tP, (c - l)aPV) > (V-, A^A^P), (A.3) 

where c > 1 was used and also the fact that aP is a projector. Since (A.3) holds for 
all quark fields in the physical subspace, this shows that the operator A^A + caP is 
bounded from below by the spectral gap of A^ A, as asserted in subsect. 5.3. 



Appendix B. Separation of bulk and boundary terms 

Equation (A.l) relates the theory in the presence of the boundaries to the one on the 
infinite lattice. The goal in the following lines is to work out the relation between the 
corresponding Dirac operators. Along the way, a separation of bulk and boundary 
terms is achieved, which allows the position-space kernels of the operators to be 
compared with each other. 
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The starting point is the identity 



R^A'^A + caP)R + (1 - R)A^ A{1 - R) = + B, (B.l) 
in which 

B = {c- l)RaPR - (1 - R)A'iAr - rA^A{1 - R) (B.2) 

denotes an operator supported at the boundaries. Since A^ A + caP operates in the 
physical subspace, and since i? is a projector, eq. (B.l) impUes 

R{A^A + caPy^^^R = R{A^A + B)-^/^R. (B.3) 



Using a well-known integral representation, the operator on the right-hand side of 
this equation can be written in the form 

(iti + S)-l/2 = 

/oo Af 
— {A^A + e)-^B{A^A + B + e)-^. (B.4) 
-oo TT 

Taken together, these equations provide a representation of [A^ A + caP)"-*^/^ (and 
thus of the Dirac operator in the presence of the boundaries) in terms of the operator 
[A'^ A)~^/'^ on the infinite lattice plus another operator localized at the boundaries. 
Both [A^ A + t"^)'^ and {A^ A + B + 1"^)'^ are in fact expected to be local operators, 
and the last term in eq. (B.4) is thus supported in the vicinity of the boundaries (up 
to exponentially small tails). 
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